To APPEAR IN The Astrophysical Journal Letters 
Preprint typeset using LMeX style emulateapj v. 26/01/00 



O 



00 
(N 

o 

u 

Oh! 

o 



> 

rn 
O 
m 



EVIDENCE FOR AGN-DRIVEN OUTFLOWS IN YOUNG RADIO QUASARS 

MiNJiN Kim' - '^ Luis C. Ho', Carol J. Lonsdale"*, Mark Lacy", Andrew W. Blain^, and Amy E. 

Kimball" 

To appear in The Astrophysical Journal Letters 
ABSTRACT 

We present near-infrared spectra of young radio quasars (Pi.4ghz « 26-27 W Hz-i) selected from the Wide- 
Field Infrared Survey Explorer. The detected objects have typical redshifts of z w L6-2.5 and bolometric lu- 
minosities ~ 10"*^ erg s"'. Based on the intensity ratios of narrow emission lines, we find that these objects are 
mainly powered by active galactic nuclei (AGNs), although star formation contribution cannot be completely ruled 
out. The host galaxies experience moderate levels of extinction. Ay « 0- L3 mag. The observed [O III] A5007 
luminosities and rest-frame 7-band magnitudes constrain the black hole masses to lie in the range ^ 10^'^- 10^^ 
Mq . From the empirical correlation between black hole mass and host galaxy mass, we infer stellar masses of 
~ 10^'-^- 10'^-^ Mq. The [O III] Une is exceptionally broad, with full width at half maximum -1300 to 2100 
km s"', significantly larger than that of ordinary distant quasars. We argue that these large line widths can be 
explained by jet-induced outflows, as predicted by theoretical models of AGN feedback. 

Subject headings: quasars: general 



L INTRODUCTION 

Supermassive black holes (BHs) and their host galaxies ap- 
pear to be closely linked in their formation and evolution (e.g., 
Magorrian et al. 1998; Gebhardt et al. 2000; Ferrarese & Mer- 
ritt 2000). While the exact origin of their relationship is still 
under debate (Kormendy & Ho 2013), theoretical models often 
invoke feedback from active galactic nuclei (AGNs) as a cru- 
cial mechanism for establishing the correlation. Outflows are 
believed to be associated with AGN feedback, which regulate 
the growth of BHs and quench star formation by expelling gas 
(e.g., Di Matteo et al. 2005). AGN feedback also appears to 
be required for cosmological models to explain the upper end 
of the luminosity function of galaxies (e.g.. Bower et al. 2006). 
Despite the potential importance of AGN feedback, direct ob- 
servational evidence of this process is limited. For example, 
extended [O III] emission in high-z radio galaxies and ultralu- 
minous infrared galaxies (ULIRGs) associated with AGNs are 
thought to be signs of AGN-driven outflows (e.g., Nesvadba 
2007;Nesvadbaet al. 2011; Harrison et al. 2012). 

Numerical simulations indicate that feedback is triggered 
during an intense and rapidly accelerated BH fueling phase dur- 
ing a gas-rich galaxy merger (Hopkins et al. 2008; Narayanan 
et al. 2010). Objects in the early feedback phase are likely to be 
both highly luminous and highly obscured by dust (Haas et al. 
2003). Thus, young quasars can be effectively searched using 
their mid-infrared (MIR) spectral energy distribution (SED). 
However, to avoid confusion with dusty starbursts, whose MIR 
SED can be similarly red, an additional diagnostic besides the 
MIR color is needed to identify AGN-powered sources. Here 
we use the presence of a strong, compact radio source to con- 
firm the AGN nature of our mid-infrared selected objects, in a 
similar fashion to Martmez-Sansigre et al. (2005; Nature, 436, 
666). Our goal is to see whether young radio quasars show sig- 



natures of outflows that might be indicative of AGN feedback 
during the early stages of their evolution. 

We use the Wide-Field Infrared Survey Explorer {WISE; 
Wright et al. 2010) to find luminous, obscured quasars (C. J. 
Lonsdale et al., in preparation). The sample is selected to be 
bright in 22 /im (^22 ^^m> ^ mJy) and to have extremely red 
colors in the MIR, 1 .25[ot(3.4 ^m)-m{A.6 ^m)] + [m(4.6 iim)- 
m(12 /im)] > 7 in the Vega magnitude system. Then we match 
the red objects from WISE to the NRAO VLA Sky Survey (Con- 
don et al. 1998) to select sources with Sia ghz/522 > 1 to 
ensure that the sample is powered by AGNs (Ibar et al. 2008). 
Finally, we only include sources with no extended radio struc- 
ture, in order to avoid evolved radio galaxies. This results 
in 156 sources. We are in the process of performing follow- 
up observations to explore the origin of these sources and to 
search for observational signatures of AGN feedback. Here we 
present near-infrared (NIR) spectroscopic data of a subset of 
these young luminous quasars. Throughout the paper we adopt 
Ho = 71 km s-i Mpc"', n„, = 0.27, and = 0.73. 

2. SPECTROSCOPIC OBSERVATIONS 

As our sources are selected from WISE and are chosen to 
be very dusty, they are expected to be very faint at shorter 
wavelengths. No reliable photometry is available in the op- 
tical or NIR. From the observed Wl (3.6 /im) flux den- 
sity and an SED template for type 2 AGNs (Polletta et al. 
2007), we estimate that our sources are fainter than 19-21 
AB mag in the K band. We obtained NIR (simultaneous J, 
H, and K) spectra for 24 objects using the Folded-port In- 
fraRed Echellette (FIRE) on the Baade 6.5 meter telescope 
at Las Campanas Observatory on 27-29 July 2012 UT. The 
data were taken in low -resolution prism mode using a 0"6 
slit. The sky conditions were clear and the seeing was 
0"8 - 1"6. The observations were taken at low airmasses 
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Table 1. The Sample 



Name 


522 urn 


'S'l.4GHz 


Exposure 


Quality 




(mjy) 


(mJy) 






(1) 


(2) 


(3) 


(4) 


(5) 


J030427.53-310838.2 


10.87 


54.99 


2536 


4 


J030629.21-335332.3 


4.83 


6.17 


1902 


4 


J035448.24-330827.5 


7.05 


7.73 


1902 


4 


J040403.61-243600.1 


4.43 


11.04 


2219 


4 


J041754.10-281654.9 


3.84 


16.32 


1902 


4 


J043921.92-315908.2 


4.57 


20.88 


2536 


4 


J130817.00-344754.3 


9.12 


68.10 


1902 


1 


J134331.37-113609.7 


3.81 


8.18 


2536 


1 


J140050.13-291924.6 


11.85 


51.92 


2536 


1 


J141243.15-202011.1 


7.41 


9.01 


2536 


3 


J143419.59-023543.6 


5.04 


36.15 


1902 


4 


J151424.12-341100.6 


7.01 


12.39 


2536 


4 


J163426.87-172139.4 


3.57 


9.49 


2536 


2 


J164107.22-054827.0 


6.26 


6.62 


2536 


4 


J165305.40-010230.5 


5.31 


12.21 


2536 


4 


J170325.05-051742.5 


11.66 


28.77 


2536 


4 


J170746.08-093916.3 


3.27 


10.86 


2536 


4 


J193622.58-335420.4 


5.27 


5.27 


2536 


4 


J195141.22-042024.5 


8.56 


20.52 


2536 


1 


J195801. 72-074609.3 


7.44 


32.79 


2536 


4 


J200048.58-280251.6 


7.19 


15.33 


2536 


1 


J202148.06-261159.2 


6.27 


6.82 


2536 


4 


J204049.51-390400.5 


4.02 


10.95 


2536 


4 


J205946.93-354134.7 


4.75 


9.13 


2536 


4 



Note. — Column (1): Source name. Column (2): Flux density 
at 22 /jm. Column (3): Radio flux density at 1.4 GHz from NVSS. 
Column (4): Integration time. Column (5): Quality of the spec- 
trum: (1) detected in [O III] A5007 and in Ha; (2) detected only 
in Ha; (3) detected emission lines are in spectral regions of strong 
telluric absorptions; (4) emission lines undetected. 

(< 1-4), with the slit aligned along the parallactic angle. Each 
object was dithered along the slit, with a position offset of 9" 
and an exposure time of 158 s in each position. The total in- 
tegration times range from 1902 s to 2536 s, split into several 
repeats of ABBA sequences. Following each science target, we 
observed a nearby AOV star (from a catalog provided by the 
Observatories) for telluric and relative flux calibration. Since a 
majority of the selected AOV stars are fairly bright, we slightly 
offset the slit from the center of the stars to avoid saturation. 
Absolute flux calibration was performed using two fainter, un- 
saturated stars that were properly centered on the slit. Since all 
of the science objects are too faint to be seen on the slit image, 
we had to apply a blind offset from the nearest Two-Micron 
All-Sky Survey (2MASS; Skrutskie et al. 2006) source in order 
to place the object on the slit. For the typical brightnesses of 
our sample, the relative astrometric uncertaintyQ between WISE 
and 2MASS is ^Q'.'4. Given the seeing conditions during the 
observing run, our absolute flux scale may have been underes- 
timated by up to 70%. Moreover, it is possible that our source 
bUnd offset acquisition procedure may have entirely missed 

*http://wise2. ipac.caltech.edu/docs/release/aIlsky/ 



some of the objects; this may partly account for the low detec- 
tion rate reported below. The FIRE spectra cover 0.8 to 2.6 yum 
with an instrumental resolution of A/ A A w 250 (1200 km s"'), 
as estimated from the full width at half maximum (FWHM) of 
the arc lines. 

We used the IDL package FIREHOSE for data reduction. 
Flat-field correction was done using images taken with an ex- 
ternal quartz lamp. We used He-nNe-nAr arc lamp spectra for 
wavelength calibration. After extracting the one-dimensional 
spectra, we corrected for telluric absorption and applied flux 
calibration using the xt ell cor package within the Spextool 
pipeline (Gushing et al. 2004; see Vacca et al. 2003 for a de- 
tailed description). 

3. SPECTRAL PROPERTIES 

Table I lists the basic properties of the sample, including the 
total exposure times and the quality of the observed spectra. 
The signal-to-noise ratio of the FIRE spectra is too low to see 
the underlying continuum. We positively detected identifiable 
emission lines in seven out of the 24 objects observed (Table 2); 
their spectra are shown in Figure 1 and 2. The redshifts lie in the 
range z « 1 .6-2.5. Six show [O III] AA4959, 5007 and Ha, but 
in J1412-2020 the lines lie in spectral regions of strong telluric 
absorption and thus are somewhat uncertain. In J 1634- 1721 
only Ha is detected. It is difficult to know why the majority 
of the sample was undetected. Possible reasons include: (I) the 
narrow-line region is heavily reddened by dust; (2) the emission 
lines coincide with strong OH lines or sky absorption; or (3) the 
slit was misplaced due to the astrometric uncertainty between 
WISE and 2MASS. In particular, since the accuracy of the tele- 
scope offsets is not precisely known throughout the night, it is 
hard to know which factors are dominantly responsible for the 
non-detection of emission lines. The detected sources are es- 
sentially indistinguishable from the undetected ones in terms of 
MIR flux, MIR color, or radio brightness. 

We do not detect broad H/3 or Ha in any of the objects. As- 
suming FWHM = 5000 km s"' for a hypothetical broad com- 
ponent and Mbh = 10^ Mq, we can use the Ha virial mass 
formalism of Greene & Ho (2005b) to estimate the flux of 
broad Ha that should be present if our objects are genuinely 
unobscured type 1 AGNs. Comparing these expected values 
with the observed upper limits of broad Ha flux, we conclude 
that the broad-line region, if present, must be reddened by 
E(B-V) > 1.2 mag. This suggests that our objects are true 
type 2 quasars rather than the lightly reddened quasars usually 
selected using NIR colors (e.g., Glikman et al. 2012). 

We follow the procedures outlined in Ho & Kim (2009) to 
fit the spectra to obtain basic parameters for the emission lines. 
We use a single Gaussian to fit H/3 and each component of the 
doublet of [O III] AA4959, 5007. The width is fixed to be the 
same for all three lines. We constrain the peak ratio of [O III] 
A4959 to [O III] A5007 to 2.88 and fix the separation between 
the two lines to the known value. For objects with no detec- 
tion of H/3, we estimate a 3a upper limit of its flux by using 
the width of [O III]. For J 1634- 1721, where [O III] is also un- 
detected, we assume a Gaussian profile with a = 500 km s"'. 
For Ha and [N II] AA6548, 6583, we simply measure the inte- 
grated flux of the lines because they are severely blended at our 
spectral resolution. 

Although H/3 is reliably detected in only two sources 
(J1400-2919 and J1951-0420), the relatively high flux ratios of 
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Fig. 1 . — FIRE spectra in the rest frame over the region sun'ounding H/3+[0 III] and Ha+[N II]. The flux density is in units of 10 " erg s ' citf -2 . The 
original spectrum and the model for the emission lines are denoted by black histogram and a thick blue line, respectively. The position of bright night sky lines is 
denoted by the shaded region. 
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Fig. 2. — Same as Fig. 1. 



[O III] A5007 to H/3 (3.4 and 7.0, respectively) suggest that the 
primary ionizing source of our sample is associated with AGNs 
rather than young stars. This finding is consistent with the 
strong radio emission relative to 22 /im flux. However, because 
we do not resolve Ha from [N II], we cannot use the [N II]/Ha 



ratio and traditional optical line ratio diagnostic diagrams (e.g.. 
Ho 2008) to completely rule out the possibility that star forma- 
tion also contributes to the ionizing continuum. For the same 
reason, it is also hard to properly constrain the extinction of the 
host galaxy from the Balmer decrement. To estimate the Ha 
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Table 2. Spectral Measurements and Derived Physical Parameters 



Name 


z 


FWHM([0 III]) 


F(H/3) 


F{\0 III]) 


F(Hq+[N 11]) 


log l^{o III] 


log -i'Ha + [NII] 


log -Pl.4GHz 


log Av 


log Mbh log Mhost 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) (12) 


J1308- 


-3447 


1.6457±0.0007 


1640±394 


<19.9 


11.3±2.3 


5.0± 3.7 


43.45±0.09 


43.13±0.29 


26.88 




8.85 - 9.61 11.3 - 12.1 


J1343- 


-1136 


2.4860±0.0009 




<10.9 


18.2±3.8 


9.8± 3.9 


43.85±0.09 


43.71±0.13 


26.34 




9.26-9.34 11.5- 11.6 


J1400- 


-2919 


1.6564±0.0006 


2120±348 


5.3±0.8: 


18.2±2.7 


47.7± 9.8 


43.66±0.06 


44.08±0.09 


26.77 


1 Q-I-O.S 
^■■^-0.6 


9.06-9.72 11.5- 12.2 


J1412- 


-2020 


1.8082±0.0056 


1900±385 


< 8.8 


11.7±9.7 


9.2± 7.3 


43.S1±0.36 


43.45±0.31 


26.09 


8.92 - 9.73 11.3- 12.1 


J1634- 


-1721 


2.0702±0.0030'' 


< 9.2 


<12.4 


25.9± 7.9 


<43.60 


43.96±0.26 


26.24 




< 9.53 < 11.9 


J1951- 


-0420 


1.5833±0.0003 


1320±211 


2.8±0.4 


19.7±2.4 


14. 5± 4.7 


43.67±0.05 


43.56±0.14 


26.33 


^■U_0.0 


9.07-9.10 11.5-11.6 


J2000- 


-2802 


2.2772±0.0028 


1870±962 


< 2.4 


12.8±5.5 


14.9±10.8 


43.66±0.19 


43.75±0.30 


26.54 


>0.2 


9.07-9.46 11.3- 11.7 



^Redshift is derived using Ha. 



Note. — Column (1): Source name. Column (2); Redshift based on [O III]. Column (3); FWHM of [O III] in units of km s^^, corrected for the instrumental 
resolution. Column (4); Flux of H/3 in units of 10^^''' erg s~^ cm~^. A measurement possibly affected by night-sky lines is flagged with a colon. Column (5); Flux 
of [O III] in units of 10"^''' erg s"^ cm'"^. Column (6): Flux of Hq+[N II] in units of 10'^''' erg s'^ cm^"^. Column (7): Luminosity of [O III] in units of erg s~^. 
Column (8): Luminosity of Hq+JN II] in units of erg s^^. Column (9): Radio luminosity at 1.4 CHz in units of W Hz^^. Column (10): Reddening (mag) in the V 
band estimated from the Balmer decrement. Column (11): BH mass in units of Mq, see text for details. Column (12): Stellar mass of host galaxies in units of iVf©, see 
text for details. 



flux from the total flux of Ha+ [N IT] , we adopt an average value 
of [N II]/Hq! (0.83 ±0.42) derived from quasars selected from 
the Sloan Digital Sky Survey (SDSS; Shen et al. 201 1). We es- 
timate Ay = 1.3 andO.Omagfor J1400-2919aiidJ1951-0420, 
respectively; the upper limit of H/3 for J2000-2802 yields a 
lower limit of Ay = 0.2 mag. These values of Ay are surpris- 
ingly low for ]VlIR-selected sources, but they might be biased 
in view of the large number of nondetections, which may be 
heavily extincted. 

Interestingly and quite surprisingly, the [O III] lines are spec- 
trally resolved for all detected sources. The FWHM of the line, 
corrected for instrumental resolution, ranges from 1320 to 2120 
km s"'. 

4. DISCUSSION 
4.1. Black Hole Mass and Host Galaxy Mass 

Since the bolometric luminosity of an AGN should not 

exceed the Eddington limit for a given BH mass (but see 
Abramowicz et al. 1988), we can estimate a lower limit to the 
BH mass from the observed luminosity, using, as commonly 
done, the [O III] luminosity as a proxy for the bolometric lumi- 
nosity. The [O III] luminosities of our sample range between 
L[o HI] ~ 10"*^ and 10"*^ erg s"'. These are likely conservative 
lower limits because we have not accounted for extinction, and, 
in any case, our observations may suffer from significant slit 
losses (see Section 2). Yet, these [O III] luminosities lie in the 
upper end of the luminosity distribution of optically selected 
z < 0.8 type 2 quasars (Reyes et al. 2008). Using a bolometric 
correction factor of L\,oi/L[o m] = 3200 (Shen et al. 201 1), we 
obtain Mbh > 10^-^- 10^-^ Mq. A caveat is that the correction 
factor can be slightly overestimated in the radio-loud AGNs due 
to the enhanced [O III] (e.g.. Stem & Laor et al. 2012). From 
the correlation between BH mass and host galaxy mass, we sus- 
pect that these systems may be hosted by very massive galaxies 
at this epoch. Taking into account the cosmic evolution of the 
BH mass-stellar mass relation (Decarli et al. 2010), we estimate 
^host > 10""'- lO"''' Mq. These mass estimates are, of course, 
highly uncertain, depending on, among other things, the strict 
vaUdity of the Eddington limit, the bolometric correction for 
[O in], reddening and slit loss corrections, and the MflH-Mhost 



relation and its uncertain redshift dependence. 

On the other hand, we can independently estimate upper lim- 
its to the mass of the hosts using the observed MIR luminos- 
ity. Given the redshift range of the sample (z« 1.6-2.5), the 
3.4 /xm band roughly corresponds to the rest-frame / band. 
Since the contribution from hot dust emission heated by the 
AGN might not be negligible in this band (e.g., Urrutia et al. 
2012), it only allows us to derive an upper limit to the stel- 
lar mass of the hosts. We convert the observed 3.4 /xm flux 
to the rest-frame / band using the template spectrum of type 
2 Seyferts from Polletta et al. (2007). Different choices in 
templates (e.g., starbursts instead of Seyfert 2s) lead to only 
moderate (^0.1 mag) differences in the resulting 7-band mag- 
nitudes. For simphcity, we adopt spectral models from Bruzual 
& Chariot (2003) by assuming a single stellar population with 
formation redshift zj = 2.5, 3, 3.5, and 5, solar metallicity, and 
a Salpeter (1955) stellar initial mass function. With these as- 
sumptions, we derive M/Lj as a function of redshift; it varies 
between 0.4 and 1. The average scatter in M/Lj at a given red- 
shift due to different Zf is 0.1 dex. Thus, Mhost < lO"-^- 
10'2-2 Mq and correspondingly Mbh < W ' - lO'-^ Mq. These 
estimates neglect the effects of dust extinction. 

The above considerations suggest that the host galaxies of 
our young quasars have stellar masses Mhost ~ lO" "'- lO'^'^ 
Mq. They populate the massive end of the galaxy stellar mass 
function at z w 2 (Daddi et al. 2004) and will evolve to become 
the most massive early-type galaxies at z = 0. 

The width of [O HI] can serve as a surrogate for the stellar ve- 
locity dispersion (cr*; e.g.. Nelson 2000; Greene & Ho 2005a; 
Ho 2009). The derived tr* range from 550 to 890 km s"' . These 
are clearly implausible: no galaxy is known to exceed tr, « 450 
km s"' (Salvianderet al. 2008). The [O HI] line must be signifi- 
cantly broadened by non-gravitational forces, presumably asso- 
ciated with an outflow or disk-wind (e.g., Whittle 1992; Greene 
& Ho 2005a). 

4.2. AGN Feedback 

There is mounting evidence that broad [O HI] emission is 
related to AGN-driven outflows (e.g., Nesvadba et al. 2007; 
Harrison et al. 2012). Figure 3 (left) shows the correlation be- 
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Fig. 3. — (left) FWHM of [O III] vs. [O III] luminosity for our sample (red circles) and various types of AGNs: high-z ULIRGs (green squares; Harrison et 
al. 2012), high-z radio galaxies (cyan triangles; Nesvadba et al. 2007, 2008), high-z type 2 quasars (blue solid stars; Nesvadba et al. 2011), high-z type 1 quasars 
(blue pentagons; Ho et al. 2012; Cano-Diaz et al. 2012), high-z compact steep-spectmm (CSS) sources {magenta crosses; Hirst et al. 2003), nearby CSS sources 
(magenta diamonds; Gelderman & Whittle 1994), and nearby type 1 quasars from SDSS (small black dots; Shen et al. 2011). The FWHMs are derived either 
from single-Gaussian (filled symbols) or double-Gaussian (open symbols) profiles. For the open symbols, the FWHMs of the broad component ai'e shown, while the 
[O III] luminosity includes both the narrow and broad components, (right) FWHM of [O III] vs. radio luminosity at 1.4 GHz. We assume fi, oc i/^-^ in order to 
estimate the radio luminosity in the rest-frame 1.4 GHz; the overall trend does not change if the sources have a flat spectrum (fu oc f"). 



tween the FWHM and luminosity of [O III] for various types 
of quasars and radio galaxies. Where the line was fit with 
two Gaussians in the literature, we plot the broad and nar- 
row components open and filled symbols, respectively; single- 
Gaussian fits are plotted with solid symbols. Clearly, our sam- 
ple has extremely broad [O III] line widths compared to other 
sources with similar redshift and luminosity. Excluding the 
SDSS quasarfl we find a weak correlation between FWHM 
and L[o lu] with a Spearman's rank correlation coefficient of 
0.54 (P = 0.0002). This correlation has been known for various 
types of AGNs (e.g.. Ho et al. 2003, and references therein), al- 
though our trend is driven by sources with extreme line widths 
(>1000 km s"') compared to those in previous studies (a few 

hundred km s"'). While the origin of the correlation between 
FWHM and L[o uij is still unclear (e.g.. Whittle 1992; Ho et al. 
2003), the very large line widths seen in our sample almost cer- 
tainly do not track gravitational motions (see §4.1). Instead, 
the most natural explanation for the supervirial velocities is 
that they reflect the kinematics of AGN-driven outflows, which 
seem to be prevalent in our sample of young radio quasars. The 
flux ratio of radio to [O III] for our sources is slightly lower 
than in high-z radio galaxies and compact steep-spectrum (CSS) 
sources, indicating that our sample is moderately radio-loud. 

A jet-induced origin for the large velocities is supported by 
the distribution of [O III] FWHM and radio power (Figure 3, 
right), which qualitatively resembles trends noted by Heckman 
et al. (1981) and Nesvadba et al. (2011). Integral-field spec- 
troscopy of high-z AGNs with broad [O III] lines finds that the 

^We include the low-z SDSS quasai's for illustration purposes only. The [O III] 
values because they were deiived from a narrow component, while the fit was done 



emission is often associated with extended outflows of ionized 
gas (e.g., Nesvadba et al. 201 1; Cano-Diaz et al. 2012). The 
energy deposition of the outflows cannot be fully explained by 
star formation, suggesting that an additional source of energy 
must come from the AGN (e.g., Harrison et al. 2012). Spatially 
resolved spectroscopy is needed to estimate the energy contri- 
bution from AGNs in our sample. 

Our W/5/i-selected, radio-emitting AGNs bears an interest- 
ing similarity to CSS sources, whose small physical extent (^ a 
few kpc) and steep radio spectrum suggest that they are young 
radio galaxies (Fanti et al. 1990). The broad and complex ve- 
locity profiles of [O III] observed in CSS sources are thought 
to arise from interaction between the radio jet and the interstel- 
lar medium of the host (Gelderman & Whittle 1994; Holt et 
al. 2008). Figure 3 shows that bona fide CSS sources tend to 
have larger [O III] widths compared to ordinary SDSS quasars; 
the line width enhancement is largest for more distant systems. 
Our sample, selected to be both dusty and compact in radio 
morphology, qualitatively resembles CSS sources, although we 
presently do not yet have sufficient radio spectral information 
to firmly confirm this association. Follow-up interferometric 
radio observations are under way. 

Theoretical models of galaxy formation that invoke AGN 
feedback often link ULIRGs and quasars through an evolution- 
ary sequence mediated by galaxy merging (e.g., Di Matteo et 
al. 2005). In this model, a massive, gas-rich merger begins 
with a ULIRG phase dominated by star formation, which then 
evolves into ULIRGs that host AGNs, young quasars, and then 
finally to optically revealed quasars with suppressed star forma- 

FWHM for the sample (Shen et al. 201 1), may be somewhat biased toward lower 
with two (naiTow+broad) components. 
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tion. The strength of AGN feedback is expected to dramatically 
increase in the later stages. In ULIRGs with AGNs, Harrison 
et al. (2012) found striking evidence for AGN feedback from 
the detection of spatially extended [O ni] emission and broad 
line widths ranging from 900 to 1400 km s"' . Their sample is 
also detected in the submillimeter, indicating substantial ongo- 
ing star formation. Interestingly, our sample has even larger line 
widths than theirs, although the stellar masses of the host galax- 



ies we deduced for our sources are comparable (e.g., Hickox et 
al. 2012). 

We thank Robert Simcoe and Rik Williams for helpful ad- 
vice on the Magellan observations. An anonymous referee 
offered valuable input. This work is supported by a KASI- 
Camegie Fellowship (MK) and the Camegie Institution for Sci- 
ence (LCH). 
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